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^ A flaw in the comparison between two diff'erent theoretical equations of state for a binary mixture 

^ . of additive hard disks and Monte Carlo results, as recently reported in C. Barrio and J. R. Solana, 

' Phys. Rev. E 63, 011201 (2001), is pointed out. It is found that both proposals, which require the 

*^ , equation of state of the single component system as input, lead to comparable accuracy but the one 

advocated by us [A. Santos, S. B. Yuste, and M. Lopez de Haro, Mol. Phys. 96, 1 (1999)] is simpler 
and complies with the exact limit in which the small disks are point particles. 
, In a recent paper. Barrio and Solana m proposed an equation of state (EOS) for a binary mixture of additive hard 
disks. Such an equation reproduces the (known) exact second and third virial coefficients of the mixture and may be 
expressed in terms of the EOS of a single component system. They also performed Monte Carlo simulations and found 
that their recipe was very accurate provided an accurate EOS for the single component system (in their case it was 
' the EOS proposed by Woodcock j|]) was taken as input. The comparison with other EOS for the mixture available 
in the literature indicated that their proposal does the best job with respect to the Monte Carlo data. Among these 
O other EOS for the binary mixture, only the one introduced by us a few years ago 1^ also shares with Barrio and 
,0, 1 Solana's EOS the fact that it may be expessed in terms of the EOS for a single component system. What we want to 
point out here is that the comparison made in Ref. is flawed by the fact that it was not performed by taking the 
^-H . same EOS for the single component system in both proposals. 

^ ' Let us consider a binary mixture of additive hard disks of diameters cti and f72- The total number density is 
' p, the mole fractions are xi and a;2 = 1 — xi, and the packing fraction is 77 = jp{a'^), where (a") = Y^i^i^i'^f ■ 
Let Z — p/pksT denote the compressibility factor, p being the pressure, T the absolute temperature, and ks the 
Boltzmann constant. Then, Barrio and Solana's EOS for a binary mixture of hard disks, Z^^{ri), may be written in 
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in 

. terms of a given EOS for a single component system, Zs{r]), as 



Z^Hv) = 1 + i(l + /Jry) (1 + C) [Zsiv) - 1] , (1) 
where ^ = (cr)^/ (cr^) and /3 is adjusted as to reproduce the exact third virial coefficient for the mixture B3, namely 



Here, 63 = 4(4/3 - VS/n) is the reduced third virial coefficient for the single component system while B3 is given by 

i 

i?3 = ^ [aiixlaf + 3ai2xlx2crt2 + ia2iXix\crl2 + a22x\(J2) , (3) 

where 



a., =Tr + 2 {allol - l) cos^i (a,/2a,,) - ^^uya} - 1(1 + l2ol)cy\ jlcyl (4) 
and Gij — (cr,; + (Tj)/2. 

The EOS for the mixture, consistent with a given EOS for a single component system that we introduced recently 
reads [01 
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We stress the fact that Eq. (||) is simpler than Eq. (|^) [which must be complemented with Eqs. In addition, 

the structure of Eq. (|^) is valid for any number of components, while the third virial coefhcient is known exactly only 
for binary mixtures. 

In Table ||, we show the results of Eqs. (|l]) and (||) when Woodcock's EOS |^ for the single component system 
[cf. Eq. (15) in Rcf. [0] is used as input in both, as well as the available MC data. Although not shown, we have 
also performed the comparison using other EOS for the single component system, namely the one by Henderson 
ours Q, and the very accurate Levin(6) approximant of Erpenbeck and Luban [Q. In all instances, it is fair to say 
that both recipes are of comparable accuracy with respect to the Monte Carlo results, their difference being generally 
smaller than the error bars of the simulation data. Nevertheless, in the particular case of using Woodcock's EOS for 
Zs{ri), as seen in Table ^, Z^^(r]) performs slightly better than Z^^(ri). This may be fortuitous since it is known 
that the Levin(6) approximant gives the most accurate approximation to Zs(j]) and when using it in Eqs. (|^) 
and the apparent (slight) superiority of Z^^{r]) is no longer there. For instance, the theoretical values of Table 
H corresponding to the packing fraction 77 = 0.6 are increased by about 0.03 when the Levin(6) approximant rather 
than Woodcock's EOS is used as input, so that in this case the accuracy of Z^^^{rj) is slightly better than that of 



TABLE I. Compressibility factor for different binary mixtures of hard disks as obtained from Monte Carlo simulations, from 
Eq. (|l|), and from Eq. (^. In the two latter. Woodcock's equation of state for the single component system is used. 
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0.20 
0.30 
0.40 
0.45 
0.50 
0.55 
0.60 



1.559(6) 
2.036(8) 
2.79(1) 
3.31(1) 
4.02(2) 
4.98(2) 
6.31(2) 



1.559 
2.040 
2.79 
3.32 
4.02 
4.97 
6.29 



1.559 
2.040 
2.79 
3.32 
4.02 
4.97 
6.28 



1.561(5) 
2.043(8) 
2.79(1) 
3.31(2) 
4.02(1) 
4.98(2) 
6.30(1) 



1.558 
2.039 
2.78 
3.32 
4.02 
4.96 
6.28 



1.558 
2.039 
2.78 
3.32 
4.02 
4.96 
6.27 



1.565(4) 
2.051(8) 
2.80(1) 
3.33(1) 
4.04(2) 
5.03(2) 
6.36(3) 



1.563 
2.048 
2.80 
3.34 
4.05 
5.00 
6.33 



1.563 
2.048 
2.80 
3.34 
4.05 
5.00 
6.33 
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0.20 
0.30 
0.40 
0.45 
0.50 
0.55 
0.60 



1.534(6) 
1.998(7) 
2.72(1) 
3.20(2) 
3.88(1) 
4.79(2) 
6.03(3) 



1.536 
1.995 
2.70 
3.21 
3.88 
4.77 
6.02 



1.536 
1.995 
2.70 
3.21 
3.87 
4.76 
6.00 



1.540(6) 
2.008(7) 
2.71(1) 
3.22(2) 
3.90(2) 
4.81(2) 
6.04(2) 



1.538 
2.000 
2.71 
3.22 
3.89 
4.80 
6.05 



1.538 
2.000 
2.71 
3.22 
3.89 
4.78 
6.03 



1.556(7) 
2.039(8) 
2.77(1) 
3.29(2) 
3.98(2) 
4.93(2) 
6.22(3) 



1.552 

2.027 

2.76 

3.29 

3.98 

4.91 

6.21 



1.552 
2.026 
2.76 
3.29 
3.98 
4.90 
6.20 



1/3 



0.20 
0.30 
0.40 
0.45 
0.50 
0.55 
0.60 



1.491(6) 
1.907(8) 
2.55(1) 
2.99(2) 
3.60(2) 
4.39(3) 



1.490 
1.905 
2.54 
3.00 
3.59 
4.39 
5.49 



1.490 
1.904 
2.54 
2.99 
3.57 
4.36 
5.44 



1.510(8) 
1.940(8) 
2.59(1) 
3.07(1) 
3.69(2) 
4.52(2) 
5.66(6) 



1.506 
1.937 
2.60 
3.07 
3.69 
4.52 
5.68 



1.506 
1.936 
2.60 
3.06 
3.68 
4.50 
5.64 



1.538(9) 
2.004(9) 
2.71(1) 
3.20(2) 
3.89(2) 
4.79(2) 
6.06(1) 



1.536 
1.996 
2.71 
3.21 
3.88 
4.78 
6.03 



1.536 
1.995 
2.70 
3.21 
3.87 
4.76 
6.00 



^Ref. § 
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Let us try to understand why both EOS give practically equivalent results. First, it may be shown that Z^^^{ri), 
while not reproducing the exact third virial coefficient B^, yields a very good estimate of it If we replace that 
estimate into Eq. (0), we get 



1 + e V 2 



By using this estimate in Barrio and Solana's EOS, we have 



where 



iUv) - 1] 



(6) 



(7) 



(8) 



According to the approximation involved in Eq. (^, the difference Z^^{r]) — Z^^{'ri) is small if the asymmetry 
of the mixture is small < 1) and/or A{ri) is small. The function A(?7) is plotted in Fig. ^ for the cases where 
Zs{ri) is given by Henderson's EOS, by Woodcock's EOS, and by the Levin(6) approximant. In all instances it is 
practically zero up to rj w 0.2 but then it grows rapidly. The most disparate mixture considered in Barrio and 
Solana's simulations corresponds to xi — 0.25, a = a^joi — 1/3, which yields ^ = 0.75. This explains the fact that 
Z^iji) - Z^^{t]) < 0.05 in the simulated cases. It should be noted however that the right-hand side of Eq. (j^) 
tends to overestimate the actual difference ^^^(77) — Z^^{ri), so that its main purpose is to illustrate the fact that 
both EOS yield practically equivalent results for not very asymmetric mixtures. On the other hand, more important 
differences can be expected for disparate mixtures, especially in the case of large densities. At a given density and a 
given diameter ratio a < 1 the smallest value of the parameter ^ corresponds to a mole fraction xi = a/{l + a) for 
the large disks, namely ^ = 4a/(l + a)^. Thus, ^ < 1 if a < 1 and, according to Eq. (R), Z^f'{ri) - Z^^^{?]) ~ A(7/). 
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FIG. 1. Plot of A(77) by assuming Henderson's EOS (dotted line). Woodcock's EOS (dashed line), and the Levin(6) approx- 
imant (solid line) for the pure fluid. 
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Let us consider now the limit in which the small disks become point particles (a 0) and occupy a negligible 
fraction of the total area. In that case, the compressibility factor of the mixture reduces to 

Zrniv)^ -^ + xiZM- (9) 

l-J] 

The first term represents the (ideal gas) partial pressure due to the point particles in the available area (i.e., the total 
area minus the area occupied by the large disks), while the second term represents the partial pressure associated 
with the large disks. In the limit a ^ with xi finite (or, more generally, xi ^ a), we have ^ xi and /3 — > 
(1 — b3/2)x2/{i + xi) [note in that limit the approximation (||) is correct], so that 

z!7''iv)'-T^+xiZsiv), (10) 



z^Hv) 



1 



1 + Xi + .X2 1 



[Zsiv) ~ 1] . 



(11) 



Therefore, while Eq. (||) is consistent with the exact property (H), Eq. (|]) violates it. In fact, the right-hand side of 
(0), with — xi, gives the deviation of Barrio and Solana's EOS from the exact compressibility factor in the special 
case a — > 0. 
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